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We consider a phase-field model where the internal energy depends on the order parameter X in 
a nonlocal way. Therefore, the resulting system consists of the energy balance equation coupled 
with a nonlinear and nonlocal ODE for X. Such system has been analyzed by several authors, 
in particular when the configuration potential is a smooth double-well function. More recently, 
in the case of a potential defined on ( — 1, 1) and singular at the endpoints, the existence of a 
finite-dimensional global attractor has been proven. Here we examine both the case of smooth 
potentials as well as the case of physically realistic (e.g., logarithmic) singular potentials. We 
prove well-posedness results and the eventual global boundedness of solutions uniformly with 
respect to the initial data. In addition, we show that the separation property holds in the case 
of singular potentials. Thanks to these results, we are able to demonstrate the existence of a 
finite-dimensional attractors in the present cases as well. 

Key words: phase-field models, smooth and singular potentials, nonlocal operators, well-posedness, 
uniform regularization properties, finite dimensional global attractors. 
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1 Introduction 

A well-known approach to study two-phase Stefan-like problems in more than one spatial dimension 
is the so-called phase-field (or difi^use interface) method. Roughly speaking, it consists in introducing 
an order parameter X whose zero level set substitutes for the sharp interface, while X = ±1 in the 
higher/lower energy phases. The classical problem is thus replaced by an order parameter dynamics, 
originated from the study of critical phenomena, coupled with the energy balance equation governing 
the temperature field. An important issue is to recover the original interface conditions and this is 
usually done by (formal) asymptotic expansions. Moreover, diffuse interface models are quite effective 
from the numerical viewpoint since there is no need for interface tracking. A significant and basic 
example of phase-field system is due to G. Caginalp (see [5], cf. also [S]), namely. 
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Abstract 



rXt AX + T4^'(X) 2-3, 



(1.1) 
(1.2) 



on a given bounded domain $7 C M'', d S {1, 2, 3}, for some time interval (0, T). Here is a rescaled 
temperature so that i? = is the equilibrium melting temperature, while £, k,T and ^ are given positive 



constants which represent the latent heat of fusion, the diffusivity, a relaxation time and a correlation 
length, respectively. The function W is (the density of) potential energy associated with the phase 
configuration. Such potential is usually a smooth double well function (typically W{r) = i(r^ — 1)^). 
However, this is just a convenient approximation of the physically more relevant logarithmic potential 
generally taking the form 

iy(r) = (1 + r) Iog(l + r) + (1 - r) log(l - r) - 7^^ 7 > 0. (1.3) 

The mathematical literature on (jl.lD - (jl.2p is rather vast and we confine ourselves to quote the pio- 
neering paper |14 and the more recent ones [T7l[T8] (see also references therein). 

In order to analyze the microscopic influences of anisotropy on the interface, in [TUl a phase- 
field system has been derived from microscopic considerations based on Statistical Mechanics. This 
system is similar to (|l.ip - (|1.2l) but for the diffusion term ^AX which is replaced by A : D^X, where 
A e W^^^ is positive definite and D^X is the Hessian of X. However, this derivation is performed 
by truncating the expansion of the interaction function (see (TO] Prop. 2.4], cf. also [11] for higher- 
order approximations). Then the author, by using formal asymptotics, deduces a modified Gibbs- 
Thompson relation in 2D. More recently, by using the same procedure, a related phase-field model 
has been obtained without approximating the interaction function (see [121 113| ). Actually, working in 
a bounded domain and choosing A = in [121 Sec. 2], the system obtained there takes the following 
form: 

+ ^Xt - kAd = 0, (1.4) 
<rVXt -K,*X + k{x)X + W'{X) = 2<;^. (1.5) 

Here K^{x) — ^"''■K{c^~^x) where <^ > is an atomic length scale and X : M'' — > R is a sufficiently 
smooth interaction kernel satisfying K{x) ~ K{—x) and such that k{x) :~ J^^ K{x — y) dy is bounded 
and nonnegative. 

The asymptotic limit \ has been analyzed in [13] and a new anisotropic interface condition 
has been obtained. On the other hand, this class of systems was already considered in some previous 
papers (cf., e.g., [U U] [5] IH] (Hj and their references). In particular, rigorous mathematical results were 
proven for smooth potentials. Well-posedness for = M and d — 1 was established in [5] through 
semigroup theory. These results were then extended to bounded d-dimensional domains with either 
homogeneous Neumann or Dirichlet boundary conditions for ?? (see [HldS])- Regarding the longtime 
behavior, the convergence of a solution to a single stationary state was shown in [15] by means of a 
suitable nonsmooth version of the Lojasiewicz-Simon inequality. Existence of an absorbing set was 
proven in 6J as well as an analysis of the w-limit sets. More recently, the results of [TS] have been 
extended to to a class of singular unbounded potentials which does not include the logarithmic ones 
(cf. [in])- Actually, in [TP equation p.5p was modified by adding an inertial term of the form aXu, 
a > 0. However, the results proved there also hold for a = 0. Then, by exploiting [19 , the existence 
of a finite-dimensional global attractor has been established in [TB]. Here we want to generalize 
such results to both smooth potentials and more general singular potentials (e.g., of the logarithmic 
type (jl.3p ). This goal is connected with the property of the solutions of getting bounded in finite 
time uniformly with respect to sufficiently general initial data. In addition, in the case of singular 
potentials, a (uniform) separation property is also needed. Here we prove all these properties for 
weak solutions originating from initial data in the energy space. In particular, in the case of smooth 
potentials, our results generalize the corresponding ones in [TS]. Moreover, for singular potentials, 
the separation property holds instantaneously, namely for t > 0, even though the initial datum is a 
pure state (see Remark 14.21 below). As a consequence, we also demonstrate the existence of a finite 
dimensional global attractor using the approach devised in [IB] . This approach exploits the only source 
of compactness for X, namely * X. Note that we cannot expect smoothing effects on X. 

For the sake of simplicity, we choose the constants in such a way that system p.4p - (|1.5p can 
be rewritten in the form 

^t + Xt- Ai") ^0, (1.6) 
Xt + J[X] + K{x)X + W'{X):^d, (1.7) 
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in n X (0, T). Here J is a linear operator which is a suitable generalization of the nonlocal convolution 
operator introduced above (see Section [2] below). Following [15], we endow the system with the 
following boundary and initial conditions 

i9 = onai7x(0,r), (1.8) 
^t=o=^Q, X\t=o^Xo inn. (1.9) 

It is worth mentioning that there are also (mainly) existence and uniqueness results for more refined 
nonlocal phase-field systems, formulated with respect to the absolute temperature, which are thermo- 
dynamically consistent also far from the equilibrium temperature (cf. [121 12S1 IMl I2S1 ED] ) ■ It would 
be interesting and challenging to extend some of the present results to such systems. 

The plan of this paper goes as follows. The main results about well-posedness and regular- 
ization properties of the solutions are stated in Section [5] Then, the corresponding theorems for 
smooth potentials and singular potentials are proven in Section [3] and Section HI respectively. The 
final Section [S] is devoted to the existence of the global attractor. 

2 Well posedness and regularization results 

In the sequel we will assume that = 1, for simplicity. We set V :— HQ{n), H L^(J7) and note 
by II • II the norm in H, by (•, •) the scalar product of H, and by || • ||p the norm in LP{n) for p G [1, oo]. 
We also let A be the Laplace operator with homogeneous Dirichlet b.c, seen either as an unbounded 
linear operator on H with domain Vr\H'^{ft) or as a bounded linear operator from V to its topological 
dual V . Then (|1.6p - (|1.8l) can be rewritten as follows: 

dt+Xt+A^^ 0, in V, (2.1) 
Xf + J[X] +/(-,X) = i9, in iJ, (2.2) 

with 

/(x,r) =/o(r)-A(a;)r, (2.3) 
where A £ L°^{n) is a given function, and 

fo e C\l;m, foir) > Vr e /, /o(0) = 0, (2.4) 
where /, the domain of /o, is an open and, possibly, bounded interval of M containing 0. We also set 

Foir) -.^ fois) ds, F(x,r) /(x,s) ds = Fo(r)-A(a;)y. (2.5) 

We assume J be a linear operator such that 

JeC{LP{n),LP{n)), \\J{u)\\Lnn) < Llluhnn), (2-6) 
for some L > independent of p and all p G [1, +oo]. Moreover we assume that 

e[l,+oo): J e CiLP'{n),L°°{n)) (2.7) 

and, finally, 

J is a compact self-adjoint operator from H to H. (2-8) 

Observe that the concrete form of the nonlocal operator J (see (|1.5D ) satisfies assumptions (j2.6p - (j2.8l) . 
provided that K is smooth enough. 

We can then introduce the energy functional 

£{i9,X):^ [ (i|,?|2 + F(-,X)-f ij[X]x). (2.9) 

It is immediate to realize that under the above assumptions £ could be unbounded from below. Thus, 
we need some condition implying that £ has some coercivity. In particular, we will consider two 
different situations. The first one deals with what we will call a smooth potential: 
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Assumption 2.1. We assume (j2.4l) with / = M. Moreover, we ask tliat 

't/kl^^' - c/ < /oWsignr < C/(|r|i+" + 1) Vr £ M (2.10) 
and for some e > 0, k/ > 0, c/ > 0, Cf > 0. 

We will speak, instead, of singular potentials in the following case: 

Assumption 2.2. We assume (|2.4p with / an open and bounded interval of M containing 0. Moreover, 
we ask that 

lim /o(r) sign r = +0O. (2-11) 

r—^dl 

Our first result deals with the case when F is smooth. In this situation, we define the energy 
space (i.e., the set of all (i?, X)'s such that £{■&, X) is finite), as the Banach space 

X := H X L^+%n). (2.12) 

Indeed, by Assumption 12. 1[ £{■&, X) is finite if and only if {■&, X) £ X. 

Theorem 2.3. Let (l2J]) - ((2^ and Assumption [Q hold and let {■&o,^o) G X. Then, there exists one 
and only one couple {■&, X) satisfying, for all T > 0, 

19 £ H\0,T;V') n L°°iO,T; H) n L'^{0,T;V), (2.13) 
X e H\0,T;H)nL°°iO,T;L'^+'{n)), (2.14) 

solving, a.e. in (0,T), system (|2.ip - (|2.2p . and enjoying the initial conditions 

^t=o = ??o, X\t=o = Xo, a.e. in ft. (2.15) 

Moreover, there exist a time Tq > depending on the "initial energy" Eq :— f (i?o, ^o) 9,nd a constant 
Co independent of Eq, such that 

Wmh^m + WmU^^n) + ||Xt(t)||L=o(o) < Co Vt > To. (2.16) 

In the case when F is singular, the energy space is given by 

X {(t9,X) £H xH : Fq{X) G L^(f7)}. (2.17) 

Actually, since the domain / of Fq is bounded, it is clear that X C H x L°°{n). Due to the constraint 
term Fq, X is not a linear space in this case. Nevertheless, it is easy to prove that has a complete 
metric structure with respect to a natural distance function (see, e.g., [281 Sec. 3] for details). 

Theorem 2.4. Let ((2^ - (|2?8)) and Assumption [2?2l hold and let {■do,Xo) G X. Then, there exists one 
and only one couple (z9, X) satisfying, for all T > 0, 

?9 e H\0,T;V')nL^{0,T;H)nL'^{0,T;V), (2.18) 
XeH\0,T;H), FoiX)eL^{0,T;L'{n)), (2.19) 
fo{X)eL^iO,T;H), (2.20) 

and solving (|2.ip - (|2.2p together with the initial conditions (|2.15p . Moreover, there exist a time Tq > 
depending on the "initial energy" Eo := £{'do,Xo) and a constant Cq independent of Eo, such that 
(|2.16p holds, together with the separation property 

<Co Vt>To. (2.21) 

A couple (i^jX) in the condition either of Theorem 12.31 or of Theorem 12.41 will be called an energy 
solution in what follows. 

Remark 2.5. An autonomous heat source term in equation (|2.ip could be easily handled (see, e.g., 
[TB]1. Some care is required in the non-autonomous case (cf. [18 for local systems). 



4 



3 Proof of Theorem [2l3 



This existence proof is a slight generahzation of 15, Thm. 1.1]. Thus we will proceed formally with 
the a priori estimates. Uniqueness goes exactly as in [15], while we will give all the details about 
(|2.16l) . In the sequel we will note with the letter c a generic positive constant, allowed to vary on 
occurrence, depending only on /o, A and L (cf. (|2.6p ). In particular, we will always assume c to be 
independent of the initial data and of time. The letter k will note the positive constants, depending 
on the same quantities as c, appearing in estimates from below. 

Energy estimate. We test (|2.ip by d, (|2.2p by Xt and take the sum. This gives (recall that J is 
self-adjoint) 

d_ 

dt 

Next, we test (j2.2p by X, to obtain 



-£{d,X) + \\Vdf + \\Xtf<Q. (3.1) 



^^ll^f + ( JM, ^) + (/(•, X), X) < X). (3.2) 
By ((2^ . Assumption \2l\ and A £ L°°(n), it is clear that 

(/(•,X),X) + (J[X],X) > K [ F{-,X) + k{J[X],X) ~ c> k\\X\\IXI + k{J[X],X) - c. (3.3) 

Jn 

Moreover, by Poincare's and Young's inequalities, 

(^, X) < i|| V^f + c||xf < i|| V^f + ^|1X||^+: + (3.4) 

for small a > and Ca depending on a. Then, summing p.ip with p.2p and taking p.3p . p.4p into 
account, we arrive at 

A(f (^?,x) + i||xf ) + n{£{^,X) + llV^f + ||X,f ) < c. (3.5) 
Integrating p.5p . we then obtain 

£{{>(t),X(t)) <c{Eoe~''' + 1), (3.6) 
for some new value of c, independent of T. In particular, by (I2.10p . this implies 

\\mr + \m\\iti<c{Eoe-^' + i). (3.?) 

Moreover, integrating p.ip over the generic interval {t,T), and using p.6p . we infer 

^(llV^f + ||X,f) <c(Eoe-''* + l). (3.8) 



Being c independent of T, the above bound can be rewritten also for T = +00. 

Regularization estimates for -d. From p.8p and the Poincare inequality, it is clear that for any 
s G [0, +00) there exists t = t(s) G [s, s + 1] such that 

||^?(r)||^ <c(Eoe-''^ + l). (3.9) 

Then, taking s > and correspondingly choosing re [s, s + 1] such that p.9p holds, testing (j2.ip by 
—Ad and integrating over (r, T) for a generic T > t, recalling (|3.5p . and using Holder's and Young's 
inequalities, we obtain 

W'^nl^ir.T-M) + I IIA^f <c(Eoe-'^- + l). (3.10) 

In particular, noting that the above holds at least for some r £ (0, 1) (corresponding to the choice 
s = 0), we have 



\\Adf <c{Eoe-^' + 1), Vl<t<T, (3.11) 
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still with c independent of T. 

Regularization estimates for X. We will now work on a generic time interval {S, S + 2) for S >\. 

Then, as a consequence of estimate (I3.10p and interpolation, we have that 

II'^IIlio((S,S+2)xO) < C. (3-12) 

Here and in what follows, C will always denote a quantity of the form 

C = g(Eoe-"^), (3.13) 

where Q is a computable nonnegative-valued monotone function, whose expression is allowed to vary 
on occurrence, depending only on the fixed parameters of the system. 

That said, we choose a sequence of small time steps t„, n G N, defined by 



27r2n2 



so 



that = - (3.14) 



and proceed by induction. Namely, we set := S and assume that, for n > 1, given > S, there 
exists tn G {tn-i,tn-i + r„) such that 

l|X(t„)||2t"l(^) <Cr-i. (3.15) 

Notice that this is surely true as n = 1 once one sets to = S, thanks to p.7p . Then, we test (12. 2p by 
|X|"'^X. In principle this would be not an admissible test function (actually at this level (12. 2p makes 
sense just as a relation in i^((0, T) x ft) and |X|"'^X needs not have the L^-summability). However, the 
procedure could be easily justified by using a truncation of |X|"'X as a test function, and then passing 
to the limit w.r.t. the truncation parameter via the monotone convergence theorem (the details are 
left to the reader). 

Moreover, since we only need a finite number of iterations, we will not take care of the 
dependence of the various constants on n and on e. We infer 

^ Mltne + (foiX), \xrx) < \xrx) + ( - J[X] + A(.)X, |Xl-X) . (3.16) 



2 + nedt 
Then, by ^M), 

Moreover, 



(/o(x),ixrx)>«:iixii^+[:+j;:-c. (3.1?) 



i^Axrx) < II |x|-x|| . ii^ii 



2 + + 



<ll^ll^ 



<fiixii^:[::jj:+cii^ii,^. (3.18) 

Finally, by dH]) and A G L°°(f7), 

( _ j[x] + A(.)x, ixi-x) < c||ixi-x||2,^iixii,±i^ < ciixii^+-^,)^ < J 11x11^+;::+;;: + c. (3.19) 

Then, integrating (|3.16p over (i„, S+2), recalling (I3.12[) . and using p. 151) and the induction hypothesis, 
we arrive at 

\\'^\\'L^(t„,S+2;L^+"-{(2)) + ll'^llL2^i"i^i'')«((t„,S+2)xn) + c||X(i„) || ^t+l(n) ^ ^"(1 + T„ ^), (3.20) 

where the second term on the left hand side ensures that condition p.l5p will be fulfilled at the level 
n + 1. In particular, the procedure can be iterated at least until n satisfies the constraint 

2 + (n + l)e „ . 8 + 9e 

^ ^ < 10, i.e., n < . 3.21 

1 + e e 



1+5 

2+(n+l)e 
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More precisely, since (8 + 9e) / e may not be an integer, our method works at least until we reach some 
nmax > 8(1 + e~^). Thus, we obtain (at least) the bound 

\W\L°°{S+l/'i,S+2-L^o+»'{Q)) + II^IIlio+9=((S+1/4,S+2)xO) < C, (3.22) 

where C additionally depends on the chosen sequence r„ . The upper bound in hypothesis (I2.10[) then 
gives also 

II/(-)^)IIl9((S+1/4,S+2)xO) <C, (3.23) 
whence, comparing terms in (|2.2p and using (|3.12l) . assumption (|2.6I) and (|3.22p - (|3.23p . we also obtain 

ll^t||L9((S+l/4,S+2)xO) < C. (3.24) 

With p.24p at disposal (indeed, any exponent p > 3, in place of 9, would be sufficient for the purpose), 
we can apply to equation (|2.ip (which is a linear PDE) a standard Alikakos-Moser iteration method 
PP with time-smoothing (see, e.g., [551 Lemma 3.5], cf. also (ISl Chap. Ill, Sec. 7]), to obtain 

II'?IIl-((S+1/2.S+2)xO) < 0, (3.25) 

for some 6 > stiU having the form 9 = Q{Eoe^'^^). 

Once we have p.25p . we can go back to (|2.2p and perform further iterations of (|3.16p . restarting 
from n — with the choice of io = -S* + 1/2. Actually, (|3.17p can still be used. On the other hand, 
thanks to (j3.25p . p.lSp can be now modified this way: 

(^,|Xrx) < |||xrx|| .+(,.+i,. ||^||^<9||X||^|-,,)^ < ^\\X\\lX<'-l'^;^ + ce'^^. (3.26) 

At this point, we can proceed with the iterations exactly as before and notice that, still in a Rnite 
number of steps, we arrive at 

ll^llL = (S+3/4,5+2;LP*(i2)) < C- (3.27) 

Consequently, using property (|2.7p . we also obtain 

ll'^MII-L~((S+3/4,S+2)xO)) < M, (3.28) 

for some constant M > having the same dependence of C. 

Remark 3.1. Unlike the case of parabolic type equations, it seems that, for equation p.2p . which has 
essentially an ODE structure, the Moser iteration method cannot be used to get directly an L°°-bound 
for X. Actually, if we try to iterate the above procedure inGnitely many times, we readily notice that 
the constants appearing on the right hand side's of the estimates (cf., e.g., p. 181) and p.l9p ) would 
explode. This is due to the fact that, while in true parabolic equations the Moser exponents grow 
exponentially with respect to n, in the present case, the growth is just linear (and, hence, too slow 
to take the constants under control). This fact also forces us to assume that the kernel J has at least 
some regularizing effect (i.e., assumption (|2.7p ). 

Once we have (|3.25p and p.28p at our disposal, we can apply a comparison principle to get the 
L°°-bound for X. Namely, we have that, for t E {S + 3/4, S + 2) and a.e. x £ ^, there holds 

Xt{t,x) + fo{X{t,x)) = X{x)X{t,x) + i}{t,x) - J[X]{t,x). (3.29) 

Let us now "freeze" the variable x. Actually, thanks to the first (|2.19p . the map t ^ X{t,x) is 
(Lipschitz) continuous for a.e. x G fl. 

Let now A > (to be chosen later) and set 

A+{x) -.^ {t e {S + 3/4,3 +2) : X{t,x)>A}. (3.30) 

Then, using ({210)) . ^2B^ and (13:29)) gives 

Xtit, x) + KfX{t, x)^+' < X{x)X{t, x) + e + M + Cf. (3.31) 



7 



for a.e. x £ fl and all t G A+(.t). Hence, dividing by X ~^'^{t, x), we obtain 

- ~X-'{t,x) + Kf < \{x)X-''{t,x) + {e + M + Cf)X-'-^+'^{t,x). (3.32) 

It is then clear that A can be taken large enough (in a way only depending on the L°°-norni of A and 
on the known constants Q, M, Cf and Kf) so that, for t G A+(x), 

-iAx-(i,x) + ^<0. (3.33) 

In particular, for such times t, the function t ^ X{t, x) is (strictly) decreasing. 

This implies that, ii t e {S + 3/4, 5 + 2) and t ^ A+(a;), then s ^ A+(j:) for aU s G [i, 5 + 2). 
In other words, if X(i, x) is smaller than A, it can never become larger than it. Thus, assuming that x 
is such that Xs(x) := X{S + 3/4, a;) > A (otherwise there is nothing to prove in view of the preceding 
discussion), integrating inequality p.33p over (5 + 3/4,t), we obtain 

X-'{t,x)>Xs{x)-' + where Xs(x) :=X(^ + 3/4, x), (3.34) 

at least for alH > S" such that X{t, x) > A. Equivalently, we can write 

At,x)<( TTT^^^^^rTT^V- (3-35) 

^ ' ^ - \KfeX%{x){t-S-3/A) + 2j ^ ^ 

Consequently, it is clear that there exists A' > 0, independent of the value o{Xs{x), such that 

X{t, x) < max{A, A'} for almost all (i, x) £ {S + I, S + 2) x n. (3.36) 

For instance, one can take A' = {SnJ^-^y/'. Of course, a similar bound from below (of the form 
X{t,x) > — max{A,A'}) can be proved in the same way. Thanks to the arbitrariness of the starting 
time S" G [1, +c») and recalling once more p.25p . we have obtained the bounds for and X in (|2.16p . 
for instance with the choice of Tq = 2 (however, see Remark 14.21 below). Then, the remaining bound 
for Xt follows by Assumption 12. 1[ (|2.6p . and a comparison of terms in (|2.2p . Theorem 12.31 is proved. 



4 Proof of Theorem 12.4 

This proof is presented in full detail (also for what concerns existence) since, to the best of our 
knowledge, this singular potential case has never been analyzed in the literature. 



4.1 Approximation 

We assume / = dom/o — (—1, 1) for simplicity. We start by approximating the singular function /□ 
by a sequence fs, S G (0, 1), such that 

fs£C\R;R), /Kr-)>OVrG/, /^(O) = 0, (4.1) 

for all 6 G (0, 1). Moreover, we require that 

Kokp-co </5(r)signr <C5(|r|2 + l) Vr G M, S £ (0,1), (4.2) 

where the constants Cs depend on 6 (and, in fact, will esplode as 5 \ 0), while kq and cq are assumed 
to be independent of S. In other words, we are asking that Assumption 12.1 1 holds with e = 1. The 
singular character of /o ensures that the lower bound in (|4.2I) can be assumed to hold uniformly in S. 
Moreover, we assume that 

/5,(r)signr < /5,(r)signr VrG/=(-l,l), VO < < (5i < 1, (4.3) 
fsir) foir) uniformly on compact sets of / = ( — 1, 1), (4.4) 
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whereas, for all r G M \ (—1, 1), fs{r) signr — > +c>o. Wc also set 

Fsir) := / fsis) ds. (4.5) 



Actually, it is easy to check that, by (14. 3p . Fg-^ < Fs^ if 62 < 5i. Moreover, we can assume that 

fs{r)r - \{x)r^ > KFs{r) - c Vr G K, (5 G (0, 1), (4.6) 

for suitable constants c > 0, k > independent of 5. The details of the construction of this ap- 
proximating family are left to the reader. For instance, one possibility could be that of taking 
fsif) ■= gs{r) + S'^iir — (1 — '5))"'")^ for r > (and an analogous choice for r < 0), where gs is 
the Yosida regularization of /o (cf., e.g., [3)- 

We also notice that, if Xq G A" (cf. (|2.17l) '). recalling that \ — 1, we have 

||>^o||p<l VpG[l,oo]. (4.7) 

Then, the above approximation permits to apply Theorem 12.31 to the system 

^s,t + Xs,t + Ads = 0, (4.8) 
Xs,t + J[Xs\ + fs{Xs) = + Kx)Xs, (4.9) 

with the initial conditions (j2.15p (note that the initial data are not approximated). This yields, for 
any 5 G (0, 1), a solution {-ds.Xs) satisfying ((2T3l) - ((2T4l) (where e = 1, cf. (|42|)), together with ((2T6l) . 
A priori, these regularity properties could depend on the approximation parameter 5. However, we 
shall see in a while that, in fact, ^-independent estimates are satisfied. 

4.2 Uniform estimates and passage to the limit 

In what follows we will assume that all constants k, c are independent of S. As we repeat the estimates 
performed in the proof of Theorem 12. 3[ it is easy to realize that, defining the approximate energy as 

£si^, X) := ^ + Fs{X) + \j[X]x) , (4.10) 

the function Sg satisfies (|3.5p with k, c independent of 6. Then, noting that 

Eo/:-f5(^9o,Xo)-^ji|??oP + i^A-(Xo)-^x2 + ij[Xo]Xo) < Eq (4.11) 

thanks to Fs < Fq, it is easy to check that the "Energy estimate" of the previous section can be 
repeated to obtain relations analogue to p.6p . p.7p . p.Sp . that hold now uniformly w.r.t. d. Moreover, 
we can test (|4.9p by fs{Xs) and use p.7p . (|3.8p . and the properties of J to infer 

II^(X^)IIl^(M+1;H) <Q(Eoe-''*), Vt>G, (4.12) 

with Q independent of S. We now show that the estimates detailed above suffice to take the limit 
5\0. Actually, (jX^ - ip^ and guarantee that, for any T > 0, 

i^s^d weakly stuT in H\0,T;V')nL^iO,T;H)nL^{0,T;V), (4.13) 
Xs~^X weakly in H^{0,T;H), (4.14) 



fsiXs)^foiX) weakly in L^(0,r;i7). (4.15) 

Here and below, we adopt the convention of overlining unidentified weak limits. Thanks to linearity 
and continuity of operator J^, it is then easy to show that, at the limit (5 \ 0, 

i^t+Xt+Ai^ = 0, (4.16) 



Xt + J[X] + foiX) - \{x)X = (4.17) 



9 



and the initial conditions (j2.15p are satisfied as well. 

Then, to conclude the proof, we have to show the identification foiX) — fo{X) almost every- 
where in (0,r) X n. To do this, we follow with some variations the argument given in [121 Sec. 2.3], 
which we report for the reader's convenience. 

First of all, letting 

W5:[0,T]^M, u;sit):^\\Xsit)~Xit)f, (4.18) 

using (|4.14p it is a standard check to verify that 

II'^<5||h1(0,T) < c. (4.19) 

Hence, we can assume that (here and below, all convergence relations are intended up to extraction 
of non-relabelled subsequences) 

Los^Lo Strongly in C°([0,T]), (4.20) 

where uj is continuous and nonnegative. Now, as a further consequence of (|4.14p . we have that 

Xs~^X in C^{[0,T];H). (4.21) 

In particular, for all t G [0,T], Xs{t) converges to X{t) weakly in H. Next, we compute the difference 
between (j4.9p and (I4.17p . test it Xg — X, and integrate with respect to the space variables. This gives 

^^11^5 - + {fs{Xs) - MX), Xs^x)^ {i^s - J[Xs] + X{-)Xs - + J[X] - A(-)X, Xs-X). (4.22) 



Now, we notice that, by (|4.12l) . (|4.14p and the first inequality in 

\\fsiXs)Xs\\L^,3^^o.T)xn)<c. (4.23) 

Consequently, 

fs{Xs)Xs -> MX)X weakly in ^^/^((O, T) x n). (4.24) 
By definition of subdifferential, we have, almost everywhere in (0,T) x Q, 

fs{Xs){Xs -X)> FsiXs) - Fs{X) > FsiXs) - FoiX). (4.25) 

Let us now test (|4.25p by a nonnegative test function G I?((0,T) x ft) and integrate. Then, by 
convexity and lower semicontinuity of the functional 

L2((0,r) X 17) ^ M, ^o(w)0, (4.26) 

J J{o,T}xn 

using (|4.15p and (|4.24l) . we obtain that 



/o(X)X0- // fo{X)X(b 
{o,T)xn JJ(o,T)xn 

= lim //" fsiXsKXs ^ X)(^ 

>liminf// (FsiXs)- Fa{X))q^>0. (4.27) 

JJ(o,T)xO 

To deduce the last inequality we have used the fact that the family of functionals {Fs}, being monotone 
increasing with respect to S going to 0, converges to Fq in the sense of Mosco (see, e.g., [3j) in 
L^((0,T) X 12). In particular, we used here the lim inf-property of Mosco-convergence: 



/ / Fo{v)(t) < liminf / / 



Fs{vs)(j> for ah vs ^ v weakly in L^((0,r) x 17). (4.28) 
T)xn 
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Thus, we have, almost everywhere in (0, T) x fi. 



/o(X)X > fo{X)X. (4.29) 
Moreover, we notice that, thanks to (j4.13p . (j4.14p . the Aubin-Lions Lemma, and assumption (|2.8p . 



{foiX) + - J[Xs] - ^ + J[XlXs - X) ^0, (4.30) 

at least in the sense of distributions over (0,T). 

Then, we can take the hmit, as (5 \ 0, of (|4.22l) . Using (|4.24l) and (I4.29p . and noting that the 
time-derivative operator is linear and continuous with respect to distributional convergence, we then 
obtain 

or, equivalently, 

~>^it)<cx^it). (4.32) 
Since w is nonnegative and w(0) = 0, we then obtain uj{t) = for all t e [0, T]. In other words, 

Xs{t)^X{t) strongly in i/, Vie[0,r]. (4.33) 
This fact, combined with (|4.14p gives 

Xs-^X (at least) strongly in i^(0,T;i7), (4.34) 



which entails in particular fa{X) = fo{X). Thus, (j4.17l) reduces to (|2.2p . as desired. 
4.3 Regularization estimates and separation property 

The above procedure is sufficient to get existence of an energy solution to our system under Assump- 
tion [521 Uniqueness of this solution is proved as in the other cases. 

To prove rigorously the separation property (I2.2ip . we go back to the <5-system (|4.8p - (|4.9p 
and start by noticing that the regularization estimates of Section [3] hold uniformly in S. Actually, 
concerning the "Regularization estimates for i?" it is easy to see that nothing changes and the analogue 
of p.9|) - (l3.11|) hold uniformly in 5. Concerning the "Regularization estimates for X", we can proceed 
as before (where we now have e = 1, of course), until we reach estimate (j3.22p . Indeed, in this part 
of the Moser iteration, we only use the estimates (|3.7l) - (|3.8p and (|3.9p - (l3.1ip . which are uniform in 
d, and the estimate from below (i.e., the first inequality) in (j4.2p . which is also independent of 6. 
Consequently, we now have the analogue of p.22p . which, for e = 1 and in the current notation, 
becomes 

ll^<5||L«=(S+l/4,S+2;Li8(n)) + \\Xs\\L^3({S+l/4,S+2)xn) < C- (4.35) 

Here and in what follows, all constants c, k and C have the same meaning as in the previous section 
and, in addition, are assumed to be independent of 6. 

However, at this point we can no longer deduce the analogue of p.23p directly, since this 
requires use of the upper bound in (j4.2p . where the constants do depend on S. 

We then have to proceed with some more care and set, for p S [1, oo), 



/ |/5(r)rsignr dr < \fs{s)\P\s\. (4.36) 
Jo 

Then, we test (14. 9p by \ f s{Xs)\^ signXg, with p to be chosen below. This gives 

^ i €i^s) + \\h{X5)\\lX\ = (|^(X^)rsignX5,^5 + \{-)Xs - J[X,]) 



Q. 



< \mXs)\ZX\+<¥^\\lX\ + \\M\lX\)- (4-37) 
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Thus, integrating in time over (r, t), where S+l/4:<T<t<S + 2, we arrive at 

- //?(^^("))+^(ll'^*ll'-^(M)xf^) + ll^*lli-^(M)xn))- (4-38) 

Let us now first choose p = 5/3. Then, according to (I4.12p . we can take r G (S" + 1/4, S + 5/16) such 
that 

\\f5iXsiTmj<C. (4.39) 

Then, by the inequality in (|4.36p . 

d^TiXsir)) < \\\fsiXs{r))\'/\JXsir)l < \\fs{Xsir))\\f\\Xs{r)l (4.40) 

and both factors on the right hand side are controlled, uniformly in S, thanks to (|4.39p and (j4.35p . 
respectively. Then, noting that the other terms on the right hand side of (j4.38p are estimated due to 
(|3?T2)) and ((435l) . we readily infer 

||/5(^<5)||LS/3((5+5/ie_5+2)xO) < C". (4.41) 

At this point, we iterate the procedure by choosing now p = 7/3 in (j4.38p . taking r G (5'+5/16, S+3/8) 
such that the analogue of (|4.39p with 8/3 in place of 2 holds (this is possible thanks to (j4.4ip ). and 
replacing (|4.40p with 

4^'iXs{r)) < \\\MXs{r))\'/\JXsir)l < (r)) Hg^g ||X,(r) (4.42) 

Thus, we finally arrive at 

\\fsiXs)\\L"i/^{{s+3/s,s+2)xn) < G. (4.43) 
A comparison of terms in (|4.9p permits now to get (in place of p.24p ^ 

ll^5,t||Li«/3((S+3/8,S+2)xn) < C* (4.44) 

and the exponent 10/3 is still sufficient to operate the Moser iteration argument with smoothing for 
•d. Thus, we arrive also in this case at the analogue of relation (I3.25p . with Q independent of S. 

With this relation at disposal, we can repeat the ODE argument of Section [3l with essen- 
tially no variation. Actually, it is sufficient to use the lower bound in (|4.2p . which is uniform in S. 
Summarizing, we have obtained 

||l9<5||L=-((l,oo)xn) < ©, ||>^<5|U°=((l,oo)xn) < M, (4.45) 

with constants Q and M independent of S. 

Remark 4.1. It is worth noting that, at the limit step, we have for free that — 1 < X < 1 almost 
everywhere (and starting from t = 0), since Fq is singular. However, (j4.45p says something more, i.e., 
that we have, for i > 1, a uniform L°°-bound independent of the approximation parameter for both 
components of the approximate solution. This is a nontrivial information especially as far as i!)s is 
concerned (indeed, we just know that the initial datum i!)q lies in H). 

As a consequence of (|4.45p . there exists (a new) A > depending only on A, J, O and Af , such that 
Xs,t{t,x) + fsiXsit,x)) < A for a.e. (t,a;) e (l,oo) X fi. (4.46) 

Moreover, 

Xs{l,x)<M a.e. in fi. (4.47) 

Let now e G (0, 1) such that /o(l - e) = 3A (cf. (|2.1ip ). Then, there exists So G (0, 1) such that 
fs{l — e) > 2A for all S G [0,(5o]- In particular, by monotonicity (cf. (|4.3p ) we have that fs{r) > 2A 
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for all r > 1 — £ and S E {0,So]. We now claim that there exists a time Tq independent of S such that, 
for (almost) all x € ^l, all t > Tq, and all S € (0, <5o], there holds that Xs{t,x) < 1 — e. Indeed, thanks 
to the above discussion, we have that, if for some S E {0,do], t > 1, and x € ^, it is Xs(t, x) > 1 — e, 
then t Xs{t,x) is decreasing. Thus, once t Xs{t,x) enters the region where Xs{t,x) < 1 — e, it 
never exits from it. Thus, freezing x as before and assuming that Xs{l, x) > 1 ^ e (otherwise there is 
nothing to prove), from (|4.46|) we have that 

Xs,tit,x) < -A, (4.48) 

at least as long as X{t, x) remains larger than 1 — e. Then, integrating (|4.46p in time and using (j4.47l) . 
we infer 

Xs{t, x) < Xs{l,x) - A{t -1) <M - A(< - 1), (4.49) 

whence 

M ~ (1- e] 

Xs{t, x) <l-e for all t > + 1 (4.50) 

and this bound is uniform with respect to 6. Proving the lower bound in the same way and passing 
to the limit S \0, we finally obtain (|2.2ip . as desired. 

Remark 4.2. Looking at the statement (and at the proof) of (|2.2ip one could think that the separa- 
tion property occurs only after some waiting time. However, the property was given in that form just 
for the sake of simplicity. Indeed, refining a bit the arguments in the proof, it is easy to demonstrate 
that (j2.21[) is in fact an instantaneous property. Namely, there holds 

- 1 + £{t) < X{t, x) < 1 - e(r) for all (t, x) e (r, +oo) x Q. and all r > 0, (4.51) 

where e(r) goes to as r \ 0. The details are left to the reader. 



5 Global attractors 

Here we establish the existence of a finite-dimensional global attractor for both smooth and singular 
potentials. The technique is the same as the one used in [TBI Proof of Thm. 4.1] for singular unbounded 
potentials. However, here we start from very general initial data in the energy space and we exploit 
the previous regularization results to define a semigroup acting on a convenient invariant set. 

Let us consider the case of smooth potentials first. From Theorem 12.31 it is clear that we 
can define a semigroup S{t) : X ^ X (cf. ^XT^) by setting [d{t),X{t)) := S'(i)(i?o, ^o), where 
("i^, X) is the unique (energy) solution to ([ST]) -([22]), ((2l^ . Note that d e C°{[Q,+oo);H) while 
X e C°([0, -t-oo); i/)nC° ([0, +oo)\ L'^^'^ [n)) . Moreover, on account of the Lipschitz continuity estimate 
[T51 (1.11)], this semigroup is also closed in the sense of [U]- Thanks to (|2.16p and p. lip . S{t) has 
an absorbing set B which is bounded in (V n L°°(f2)) x L°°{Vt). Hence we can find to > such that 
S{t)B C B for all t > t^. Then, without loss of generality, we can suppose that Iq — Q and assume 
that B is an invariant set for S{t). Moreover, we can endow B with the V x i?-metric and obtain a 
complete metric space X. 

Wc can then prove the following 

Theorem 5.1. Let the assumptions of Theorem 12.31 jioJd and suppose that 

Ao := ess sup A(a;) < 0. (5.1) 

Then the dynamical system (X, S{t)) has a finite-dimensional connected global attractor. 
For readers' convenience, we report here below the argument of [16j . 
Proof. Consider (??oi,^oi) £ X, i ^ 1,2, set 

{m,m) = {{^i-^2Kt),{x,-x2){t)) 
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where {{)i{t),Xi{t)) = S{t){doi,Xoi) for t > 0, and observe that 

dt+Xt + Ad^O, a.e. in (0, +00) X 17, (5.2) 
+ J[X] + /(•, Xi) - /(•, X2) = a.e. m (0, +00) x n. (5.3) 

Let us multiply equation (15. 2[) by ^z?(t). Integrating over we get 

from which, using the Young and Poincare inequalities, we derive the estimate 

^llV^f + «i||Vi?f <||x,|p. 
and, by comparison in (j5.3l) . we deduce 

A||Vz?f + «i||V^f < ci (||xf + II J[X]f + ll^^lp) , (5.4) 
for all t > 0. On the other hand, multiplying (j5.3l) by X{t) and integrating over we find 

^^ll^f + - f{;X2),X) = -{J[X],X) + (^,X). 

Thus, on account of ([O)) . ((^ and ((??T|) . we get 

1 d 



|^l|X|r + Ao||X|r <-(J[X],X) + (^,X). 



Then, Young's inequality gives 

~\\xr + ^\\xr < c,A\\j[m' + imn ■ 

Since J is compact and self-adjoint on H, we can find a finite-dimensional projector such that 

\\j[v]r<^j\vr+c\\u,Mr, (5.5) 

for all V G H. As a consequence we have 

^^llXf + ^||Xf <c.„(||n[X]|p + ||,9||^). (5.6) 

Adding inequality (|5.4p multiplied hy fi — to (j5.6p and using the analogue of (j5.5p to estimate 
the term || J[X]|p on the right hand side of (j5.4l) . we infer 



d_ 



(5.7) 



Vliv^f + i||xf ) + ^iMllv^f + y l|x|p < (||n[x]f + . 

Therefore, from (j5.7p . we deduce 

||^(i)fy + ||X(i)f < c,„e-'=-o* (||^(0)||2. + ||x(0)|p) (5.8) 
+ CAo r (ll^(r)f + ||n,JX(r)]|p)dr, 



1/2 







for all t e [0, T] and any fixed T > 0. 

We now introduce the following pseudometric in X 



T 

dT((^oi,Xoi),(i?o2,Xo2)) = ( / (11(^1 -^?2)(r)f + ||nAo[(Xi-X2)(r)]f)dT 
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and we recall that a pseudometric is (pre)compact in X (with respect to the topology induced by the 
X-metric) if any bounded sequence in X contains a Cauchy subsequence with respect to (see, for 
instance, i21j Def. 2.57]). 

Let us prove that dx is precompact in X. Let {(i^onj^On)} C X (X is bounded) and set 
(■!?«(<): ^™(i)) = S{t){don,Xon)- Thauks to we have that {??„} is bounded in L'^{Q,T;D{A)) n 

H^{0,T; H). Therefore it contains a subsequence which strongly converges in L^{0,T,V). On the 
other hand, we have that [X„]} is bounded in L°°{0,T; H). Also, by comparison in (j2.2p . we 
deduce that {(X„)t} is bounded in L°°(0, T; H). Therefore {{Ux^ [X„])t} is bounded in L°°(0, T; H) as 
well. Then {nAp[X„(-)]} contains a subsequence which strongly converges in L^{0,T; H). Summing 
up {(t^onj^ori)} contains a Cauchy subsequence with respect to d^. 

From (j5.8p . we deduce that there exists t* > such that 

\\S{n{^01,Xoi) - S{t*){^02,Xo2)\\x 

< - ^02,Xoi - Xo2)\\x + c^„dt,{{'^ouXol),{^02,Xo2))■ 

ilence S{t) has a (connected) global attractor (see [H] Thm. 2.56, Prop. 2.59]) of finite fractal (i.e., 
box counting) dimension (cfr. [20, Thm. 2.8.1]). I 

In the case of singular potentials, we can first notice that, by (j2.20l) and the first (j2.19D . t H> /^^ ^o('*'^(^)) 
is absolutely continuous in [0,T] due to [Tl Lemma 3.3, p. 73]. Hence, taking X defined by (|2.17p as 
phase space, the semigroup S{t) defined as above takes X to itself for all t > 0. Using again (j2.16p and 
on account of the separation property (|2.2ip . it is not difficult to realize that there exists an absorbing 
set of the following form (note that (13. lip still holds): 

B{R,/3) := {{^,^) : ]]i9|ly < R, -I + /3 < ip < I - ^ a.e. in n}, (5.9) 

for a suitable pair of constants {R, f3) e (0, +oo) x (0, 1). Then, reasoning as above, we can suppose 
that B{R, (3) is invariant for S{t) and we can endow it with the V x iJ-metric. The resulting complete 
metric space X is now our phase space and we have 

Theorem 5.2. Let the assumptions of Theorem 12.41 and (|5.ip hold. Then, the dynamical system 
{X, S{t)) has a Enite dimensional connected global attractor. 

The proof goes as above. 

Remark 5.3. Assumption (|5.ip is crucial for investigating the asymptotic behavior (see, e.g., ]W[ 
Thm. 1.2, (1.19) and Rem. 3], Thm. 4.6 and (A4)] and [H Ass. 4, Sec. 2 and Lemma 3.1]). Also, 
we recall that if /o is real analytic then, following [151 [H] , one can prove that the w-limit set of any 
pair (i!?o,Xo) in the energy space X reduces to a singleton {(0,Xoo)}, where Xq solves the stationary 
problem 

-'^[^cx,] + fo{Xoo) - X{x)Xoo = 0, a.e. in n. 
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